In the Title 

Kindly replace the Title with the following: 
A-METHODS OF ESTIMATING PRECIPITATION CHARACTERISTICS 



2 



In the Specification 
Kindly replace paragraphs [0001] through [0021] with the following: 

Related Application 

This is a §371 of International Application No. PCT/FR2004/002692, with an international 
filing date of October 21. 2004 (WO 2005/040853 Al, published May 6, 2005), which is based on 
French Patent Application No. 03/12299, filed October 21, 2003. 
Technical Field 

Th e pr e s e nt This invention relates to armethods of estimating precipitation characteristics and A 
in one p articula r aspect, precipitation rate for solid precipitation. 
Background 

It is known from th e stat e of th e art that the characteristics of a radar image can be used to 
estimate rainfall characteristics. In particular, Europ e an Pat e nt EP 1 049 944 describes a technique 
for estimating rainfall using a radar. In that patent, provision is made for the following steps: 

[[-]] using thata dual-polarized radar to measure the differential phase shift (Odp) and the 
apparent reflectivity Z in at least one of the H (horizontal) or V (vertical) polarizations and 
over a given range £{ro, ri}) of path radii relative to saidthe radar; 

[[-]] determining an estimate of the value No* representative of the size distribution of the 
rain drops, on the basis of the differential phase shift difference over the range ro to ri, and on 
the basis of an integral of a function of the apparent reflectivity Z, over the range £{ro, ri}); 
and 

[[-]] deducing the value of the precipitation rate at a point on the basis of No* and of the 
apparent reflectivity at saidthe point. 
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A method of estimating rainfall is also known from PCT Pat e nt WO_03/007016^ which de- 
scribes a method of estimating a precipitation rate by means of a dual-polarized radar, that method 
being characterized by the following steps: 

[[-]] using thata dual-polarized radar to measure the differential phase shift Odp and the 

attenuated reflectivity Z in at least one of the H or V polarizations and over a given range 

f(rl , rO}) of path radii r relative to saidthe radar; 

[[-]] determining an estimate of the value K(rO) of the specific attenuation at rO on the 
basis of the profile of the attenuated reflectivity measured in this way, and on the basis of the 
differential phase shift difference over the range rO to rl ; determining an estimate K(r) of the 
specific attenuation at r as a function of the attenuation K(rO) determined in this way and of 
the attenuated reflectivity profile Z(r); and determining the precipitation rate R(r) once K(r) is 
known. 

Those various solutions make it possible to characterize liquid precipitation and te-estimate 
precisely rainfall rate (in millimeters per hour (mm/h), but they do not make it possible to estimate 
the characteristics of solid precipitation such as snow. 
Summary 

This invention relates to a method of estimating precipitation characteristics including ac- 
quiring a radar image including at least a vertical plane of a precipitation zone, processing a vertical 
profile to generate digital signals representative of reflectivity in vertical direction z, integrating the 
signals representative of reflectivity by assimilation of a reflectivity vertical profile in an aggregation 
model to generate a signal representative of the profile in the vertical plane of a mean particle diam- 
eter weighted by mass of each particle, and determining concentration of the solid particles on the 
basis of signals previously determined. 
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An object of th e inv e ntion is to provid e a nov e l m e thod that r e m e di e s that drawback. 

To this e nd, and in its most g e n e ral acc e ptation, the inv e ntion provid e s a m e thod of esti - 
mating pr e cipitation characteristics and in particular pr e cipitation rat e for solid precipitation, said 
method comprising an acquisition and proc e ssing st e p consisting in acquiring a radar imag e includ 
ing at l e ast a vertical plan e of a pr e cipitation zon e , and in proc e ssing a v e rtical profile so as to d e liv e r 
digital signals r e pr e s e ntativ e of r e flectivity in th e v e rtical direction h, said m e thod boing char 
acteriz e d in that it furth e r compris e s an integration st e p consisting in integrating said signals rep - 
r e s e ntativ e of r e fl e ctivity so as to d e liv e r a signal r e pres e ntative of th e profil e in th e v e rtical plan e of 
a th e m e an particl e diam e t e r w e ight e d by th e mass of e ach particl e , and a d e termination st e p con 
sisting in d e t e rmining th e conc e ntration of th e solid particles on th e basis of th e signals computed in 
th e pr e ceding st e ps. 

Pr e f e rably, th e int e gration step consists in d e t e rmining th e variable Z(h) of th e radar obs e rv 
abl e in mm /m as a function of th e altitud e h on th e basis of said radar imag e , and in d e t e rmining 
said m e an diameter Dm(h) of th e particl e s by r e solving the following e quation: 

fD m b-5 fll/Z ^ 

- 0.25k egaB— W^Z+-J l© m (2) 

fh m \6Zfh J 

wh e r e : 

Z is th e radar obs e rvabl e to b e inv e rt e d in mmV; 
B m is in m e t e rs (m); 

a and b ar e coeffici e nts sp e cific to particles of th e "aggregat e " type; for e xampl e , th e co e ffi 
cient a is e qual to 3518 4 and th e co e fficient b is e qual to 3.16; 

keff is th e co e ffici e nt of e ff e ctiv e n e ss of the aggr e gation proc e ss to b e adjust e d, said coeffi 
cient ke g b e ing e qual to 0.3. 
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Int e gration of (2) r e quires an int e gration boundary condition. Advantag e ously, said inte - 
gration boundary condition is d e termin e d so that th e valu e D m (h) at th e top of the cloud corr e sponds 
to th e predet e rmin e d valu e for the total numb e r of particl e s at th e top of tho cloud. 

In an advantag e ous impl e m e ntation, th e profil e of th e total numb e r of particl e s 114(h) is d e ter 

min e d by th e following equation: 

ft,(h)-x.Z(h)/D ffl {h) 6 
i g 

wh e r e x is e qual to 25. 4 10 - 

In another implem e ntation, th e m e t e orological param e t e r No(h) is d e t e rmin e d by th e follow 
ing e quation: 

^o(h) y.Z(h)/D m {fe) 
wh e r e y is equal to 102 ICT ^t 

In a third impl e m e ntation, th e meteorological param e t e r corresponding to th e profile of the 
ic e wat e r cont e nt IWC(l%) (in g/m * ] is d e t e rmin e d by th e following e quation: 

IWC(h)-wZ(h);0) ffl (h) 3 
wh e r e w is equal to 1.25 lO" 4 ^ 

In a fourth implem e ntation, the m e teorological param e t e r corresponding to the profil e of th e 
solid pr e cipitation rate R(h) (mm/h e quival e nt m e lt e d) is det e rmin e d by th e following e quation: 

R(h)-r.Z(h>0) ffl (h)^ 
wh e r e r is equal to 4.698 10 ^°t 
Brief Description of the Drawings 

: Fh eSelected aspects of the invention will be better understood on reading the following de- 
scription, given with reference to a non-limiting impl e mentation. example and selected graphic rep- 
resentations below: 
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Figur e 1 shows an e xampl e Fig. 1 is a graph of a vertical profile of Z to be inverted (in this 
example, the isotherm 0°C is at ground level); 

Figur e 2 shows th e Fig. 2 is a graph of a profile of Dm resulting from Z being inverted by 
the aggregation model compared with thea conventional estimator; 

Figur e 3 shows th e Fig. 3 is a graph of profiles of NO and of nT resulting from Z being in- 
verted by the aggregation model compared with thea conventional hypothesis and obser- 
vations; 

Figur e 4 shows Fig. 4 is a graph of the vertical profile of Z to be inverted; 

Figur e 5 shows Fig. 5 is a graph of the profile of Dm resulting from Z being inverted by the 
aggregation model compared with the conventional estimator; 

Figure 6 shows Fig. 6 is a graph of the profiles of NO and nT resulting from the Z being 
inverted by the aggregation model compared with thea conventional hypothesis and obser- 
vations; 

Figure 7 shows Fig. 7 is a graph of the profile of R resulting from Z being inverted by the 
aggregation model compared with thea conventional estimator; 

■ Figur e 8 shows Fig. 8 is a graph of the sensitivity of the retrieval of the equivalent precip- 
itation rate (mm/h) to the particle density law p(D) oc D-y; and 

Figur e 9 shows Fig. 9 is a graph of the sensitivity of the retrieval of the vertical Doppler 
velocity to the particle density law p(D) oc D-y. 
Detailed Description 

This disclosures provides methods of estimating precipitation characteristics and, in a par- 
ticular selected aspect, a precipitation rate for solid precipitation, the method comprising an acqui- 
sition and processing step comprising acquiring a radar image including at least a vertical plane of a 
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precipitation zone, and in processing a vertical profile to deliver digital signals representative of 
reflectivity in the vertical direction h, an integration step comprising integrating signals repre- 
sentative of reflectivity to deliver a signal representative of the profile in the vertical plane of a mean 
particle diameter weighted by the mass of each particle, and a determination step comprising deter- 
mining the concentration of the solid particles on the basis of the signals computed in the preceding 
steps. 

Reference will be made to the drawings below as appropriate. Steps in accordance with se- 
lected aspects of the method may be seen in one sense by reference to Figs. 1-3. 

In a preferred aspect, the integration step comprises determining the variable Z(ti) of the radar 
observable in mm 6 /m 3 as a function of the altitude h on the basis of the radar image, and determining 
the mean diameter Dm(h) of the particles by solving the following equation: 

/Dm b-5 fll/Z ^ 

= -0.25keffaD 10' 18 Z+ I |D m (2} 

fh ™ l6Zfti J 

where: 

Z is the radar observable to be inverted in mm 6 m" 3 ; 
Dm is in meters (m); 

a and b are coefficients specific to particles of the "aggregate" type; for example, the coeffi- 
cient a is equal to 35184 and the coefficient b is equal to 3.16; 

fc ff is the coefficient of effectiveness of the aggregation process to be adjusted, said coeffi- 
cient kg f r being equal to 0.3. 

Integration of (2) uses an integration boundary condition. Advantageously, the integration 
boundary condition is determined so that the value DUh) at the top of the cloud corresponds to the 
predetermined value for the total number of particles at the top of the cloud. 
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In an advantageous implementation, the profile of the total number of particles rit(h) is deter- 
mined by the following equation: 

n x (h) = x.Z(hyD^ h) 6 
where x is equal to 25.4 10~ 18 . 

In another implementation, the meteorological parameter Nn(h) is determined by the follow- 
ing equation: 

NnOi) = vmVDn fli) 7 
where y is equal to 102 10~ 18 . 

In another aspect, the meteorological parameter corresponding to the profile of the ice water 
content IWCfhK'm g/m 3 ) is determined by the following equation: 

IWOh) = wZrhVD^ 3 

1 9 

where w is equal to 1.25 10" . 

In yet another aspect, the meteorological parameter corresponding to the profile of the solid 
precipitation rate R(h) (mm/h equivalent melted) is determined by the following equation: 

R(h) = r.Z(h)/D m fh) 235 
wherein r is equal to 4.698 10~ 10 . 

The method applies e ss e ntially to stratiform precipitation. It considers that the-ice-forming 
nuclei are activated only at highly negative temperatures, i.e.,, at the top of the cloud. Th e ic elce 
crystals formed at high altitude settle and grow as they fall, either by sublimation of the ambient sat- 
urating water vapor, or by collection and freezing of supermelted cloud water droplets, or eke-by 
random aggregation of their collisions with other ice crystals. Of the three growth processes, only 
aggregation changes the ice particle concentration. The characterization method is based e ss e ntially 
on a simplified description of the aggregation mechanism. Th e inv e ntion is bas e d on We use a class 
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of "profiler" methods in-that ^inverts the vertical profile of reflectivity measured in the ice , so as to 
extract therefrom the vertical profile of the solid precipitation rate. 

TheSelected steps of th e inv e rsion m e thod ar e as follows include : 

1 - Particle size distribution expressed in "equivalent melted diameter" is assumed to be 
exponential, i.e.: 

N(D)=N 0 exp(-4D/D m ) (1) 
where N(D) is the concentration of particles per cubic meter (m 3 ) and per diameter range, and 
No and D m are the two parameters that characterize distribution. 

2 - The top h max and the base h min of the layer of solid precipitation are determined; 

a. h max is the maximum altitude of the measured reflectivity profile Z(h). 

b. hmin is either the altitude of the isotherm 0°C if the ground-level temperature is 
positive, or it is ground level if the ground-level temperature is negative. 

3 - The profile of the parameter Dm(h) in the range h max to h min is then determined by resolv- 
ing the following differential equation: 

/Dm b-5 fll/Z ^ 

= -0.25ketfiD 10 -18 Z+ I |D m (2) 

fh m l6Z/h J 

where: 

- Z is the radar observable to be inverted inmmV 3 ; 

- D m is in meters (m); 

- a and b are coefficients specific to particles of the "aggregate" type, equal respec- 
tively to 35 1 84 and to 3. 1 6 on the basis of the observations of J.D. Locatelli and P.V. 
Hobbs , Fall speeds and masses of solid precipitation particles. J. Geophvs. Res., 79, 
2185 - 2197 (1974) , the subject matter of which is incorporated herein by reference ; 
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- keffis the coefficient of effectiveness of the aggregation process to be adjusted (the 
value k e ff= 0.3 seems correct). 

4 - The integration of (2) takes place from the top, where the boundary condition is ex- 
pressed by fixing the total number of particles ht {or the number of ice-forming nuclei acti- 
vated at the top of the cloud). It is possible to take n T (h max ) = 10 6 m~ 3 , which makes it pos- 
sible to express the boundary condition D m (h maj( ) as: 

Dm(H max )=25.4 W' I8 (Z(h max /n 7 {h m J) ,/6 (3) 

5 - Once the profile D m (h) from h max to h m j n has been determined, the profiles of the other 
parameters of interest are computed by the following expressions: 

a. Profile of TV* N 0 (h)=102.10 12 Z(h)/D m (h) 7 

b. Profile of the total number of particles n T (h) ftm m" 3 }): N T (h)=102.10" 12 Z(h) 
/D m (h) 6 

c. Profile of the ice water content IWC(h) (in g/m 3 }): 

IWC(h)=l .25. 10" 12 Z(h)/D m (h)3 

d. Profile of the solid precipitation rate R(h) (mm/h equivalent melted). By using the 
terminal fall velocity determined by Locatelli and Hobbs for aggregates: {{vt=107.6 
D° 65 (D in m)}), R(h) is expressed by: 

R(h)=4.698. 10- 10 Z(h)/D m (h) 2 35 

The following description relates to another versie nselected aspect of tfeea method of pro- 
cessing for determining the precipitation rate. 

For rainfall, the algorithm used is the ZPHI algorithm which is the subject of French Pat e nt 
Docum e nts FR 98/00714 and FR 01/09206 , the subject matter of which is incorporated by reference . 
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For solid precipitation, the nov e l algorithm that has just b ee n d e v e loped is the subj e ct of tho 
following descriptio n applies . 

This algorithm The method for estimating the precipitation rate for solid precipitation b e longs 
te4heincludesa class of "profiler" algorithms ia-that it-inverts the vertical profile of reflectivity mea- 
sured in the ice , so as to extract therefrom the vertical profile of the solid precipitation rate. 

The algorith mm ethod applies e ss e ntially to stratiform precipitation. It considers that the ice- 
forming nuclei are activated only at highly negative temperatures, i.e.^ at the top of the cloud. The 
ice crystals formed at high altitude settle and grow as they fall, either by sublimation of the ambient 
saturating water vapor, or by collection and freezing of supermelted cloud water droplets, or else by 
random aggregation of their collisions with other ice crystals. Of the three growth processes, only 
aggregation changes the ice particle concentration. The algorithm is based essentially on a simplified 
description of the aggregation mechanism. TheSelected steps of th e inv e rsion m e thod are as follows: 
Particle size distribution expressed in "equivalent melted diameter" is assumed to be 
exponential, i.e.: 

N(D)=N 0 exp(-4D/D rn ) (1) 
where D is the equivalent melted diameter of the ice particle; 
N(D) is the concentration of particles per m 3 and per diameter range; and 
N 0 and are the two parameters that characterize distribution. 

The top h m ax and the base h m j n of the layer of solid precipitation are determined; 
h max is the maximum altitude of the measured reflectivity profile Z(h). 
h min is either the altitude of the isotherm 0°C if the ground-level temperature is 
positive, or it is ground level if the ground-level temperature is negative. 
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The profile of the parameter D m (h) in the range h max to h m j n is then determined by resolving 
the following differential equation, whose solution can be determined analytically: 
dD m 4 d keffJ(b,d)a N «) 1 dZ 



= -24 ZD + D m (2) 

5h T(4+d) a(a-l) m (a-l)Zdh 

where: 

foo fx r i 

J(b,d) = | exp(-4x)x b | I x d -x ,d |exp(-4x')dx'dx 

Jo Jo L 

describes the collision frequency; 

Z is the radar observable to be inverted inmm 6 m- 3 : 
D m is in m; 

a, b t c, and d are coefficients dependent on the density law F(D) for ice particles, which 
law is assumed to vary at D y (where Yean take values in the range 0.25 to 1 . 1 depending 
on the type of particles). These coefficients are drawn from Mitchell's Theory (Journal 
of Atmospheric Sciences, 53, 12, 1996) , the subject matter of which is incorporated by 
reference, for representing the following power laws: 

v t (D) = cD d (where v t is the terminal fall velocity of the particle of equivalent melted 

diameter D): 

A = aD b (where A is the effective cross-section of the particle of equivalent 
melted diameter D); 

With coefficients depending on frequency and anas a function of (D), representing by 
a power law the relationship between reflectivity (measured parameter), the concentration of 
the particles, and their mean equivalent diameter: 

Z = aN 0 D m d 
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kefr is the coefficient of effectiveness of the aggregation process to be adjusted (the 
value kefr= 0.3 seems correct immediately above the isotherm 0°C). 
Kindly replace paragraphs [0023] through [0026] with the following: 
Once the profile D m (h) from h ma x to h m ; n has been determined, the profiles of the other para- 
meters of interest are computed by the following expressions: 
Profile of N 0 : 
Z(h) 

N 0 (h)= 

dD m (hf 

Profile of the total number of particles n T (h) £{in m" 3 }): nj{h)=Q.25No(h)D m (h) 
Profile of the ice water content IWC(h) (in g/m 3 }): 

IWC(h)=(l.2272)l0 4 N 0 (h)D m (hf 
Profile of the solid precipitation rate R(h) (mm/h equivalent melted): 

T(4+d) 

J?ftHl-885)10 6 c N 0 (h)D m (h) (4+d) 

4 (4+d) 

Vertical profile of the Doppler velocity V D (h) (in meters per second (m/s)), given by a 
power law dependent on density _(D): 
V D (h)=pD m (hf 

Because of the dependency of the coefficients relative to the density law F(D) 9 the retrieval of 
the precipitation rate depends critically on the parameter 7, as shown by Figur eFig. 7. The density 
law thus constitutes a key parameter. 

By operating in vertical firing, it is possible to measure the Doppler velocity Vo(h), and te 
compare it with Vo(h) computed by the algorithm. The Doppler velocity profile is discriminating, as 
shown by FigweFig. 8. By successive adjustments, it is thus possible to determine the parameter y 
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of the law F(D) for which consistency between V D (h) as computed and V D (h) as measured is 
achieved. FigufeFig. 4 shows an example of a profile inverted by the algorithm (in thisthat example, 
the isotherm 0°C is at ground level). FigweFig. 5 shows the<:urves resulting from application of the 
aggregation model and of th ea conventional model, making it possible for the profile of D M resulting 
from the inversion of Z by the aggregation model to be compared with the conventional estimator. 

FigufeFig. 6 shows the-curves of the profiles as a function of altitude (m) above the isotherm 
0°C. It enables the profiles of No and n T resulting from Z being inverted by the aggregation model to 
be compared with the conventional hypothesis and observations. FigweFig. 7 shows the profile of R 
resulting from Z being inverted by the aggregation model compared with the conventional estimator. 
Figur e Fig. 8 shows the sensitivity of the retrieval of the equivalent participation rate (mm/h) to the 
particle density law p(D) oc D' Y . FigureFig. 9 shows the sensitivity of the retrieval of the vertical 
Doppler velocity to the particle density law p(D) oc D' Y . 
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